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£T) ' Abstract. Two decades ago Matsui and Satz suggested that Debye screening in the 

■ quark-gluon plasma would result in J/ip suppression in heavy ion collisions. Much has 

happened in the subsequent years, and the picture of quark-gluon plasma at present 
is rapidly evolving - what does it imply for the J/ip suppression? What are the recent 
s * i RHIC and SPS results trying to tell us? What else has to be done? This talk is an 

| attempt to address these questions. 
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1. Introduction 

This week thirty two years ago the physics world was shaken by what has later become 
known as "November revolution": the discovery of a new narrow heavy resonance at 
BNL and SLAC [TJ [2]. This revolution (unlike many others) is guaranteed to have a 
lasting effect - first, it added a new "charmed" quark to the list of fundamental building 
blocks of Nature; second, the explanation of the narrow width of the new resonance 
was one of the first successes of QCD - the newly born theory of strong interactions. 
The large mass of the charmed quark m c ~ 1.3 GeV and the smallness of the running 
QCD coupling constant at this scale a s {m c ) ~ 0.3 allowed to build a successful theory of 
"charmonium" decays. Moreover, the size of charmonium R ~ (m c <C Aqcd 

is sufficiently small to justify an attempt of using a perturbative Coulomb potential to 
describe the bound state structure. 

All of this makes charmonium, and heavy quarkonia in general, an extremely 
attractive starting point for an investigation into the structure of hadrons. By studying 
the distortions of the naive " QCD hydrogen atom" picture caused by non-perturbative 
effects we glean information about the structure of the vacuum. At high temperature, 
following Matsui and Satz |3], we hope to be able to understand the nature of quark- 
gluon plasma by looking at modifications of the spectral density of heavy quark 
correlation functions; experimentally, we have an access to them by measuring the 
nuclear modification of charmonium and bottomonium production. 
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2. J/ip in QCD vacuum and the origin of confinement 

Confinement effects in heavy quarkonia are often described by adding a linear term to 
the Coulomb potential, 

V(r) = -^ + ar, (1) 
6 r 

where a is the string tension; the resulting "Cornell potential" [I] has been shown 
to describe the spectrum of charmonium states quite well (for a recent review, see 
[5]). Another strategy of accounting for non-perturbative effects is to use the operator 
product expansion (OPE) and consider the leading power correction to the correlation 
function of the heavy quark-anti-quark current [6] stemming from the dimension four 
"gluon condensate" (a s G 2 ). This procedure is analogous to considering a correction to 
the energy levels of a hydrogen atom in an external electric field E 2 . Because of the 
color neutrality of quarkonium, the leading contribution arises from quadratic Stark 
effect [3 |Sj, and the corresponding shift of the energy levels SE ~ (a s r 2 E 2 /e B ), where 
€b is a typical binding energy. In QCD, the resulting energy shift is positive - it therefore 
corresponds to a confining interaction but a one quadratic, not linear, in the inter-quark 
distance. 

This problem can be seen in the formal OPE expansion for the effective heavy quark 
potential [9] (see jTUl (TT] for reviews) : 

^V(r) c - {N l~ c 1)a "l r) (l +E«n<(r) + c 3 A^ CD r 3 ) , (2) 

where the second term in the brackets describes the higher order perturbative 
corrections, and the third one is the leading power correction giving rise to the quadratic 
term in the potential. Since the linear behavior of confining potential has been firmly 
established by now in lattice calculations, we have to understand the origin of the linear 
confinement from the point of view of OPE. This is necessary because OPE provides the 
link to partonic description, and we need it for addressing the problem of quarkonium 
interactions in matter. Several interesting scenarios have been proposed to achieve this 
goal, and space does not allow me to describe them in this talk (see e.g. [TUl, [TTJ ) ; all of 
them however involve a modification of the gluon propagator, or the QCD coupling, in 
the infrared region (the necessity of such a modification is made evident by the Gribov 
copies problem [12] - the existence of multiple solutions of gauge fixing condition at 
large coupling). In terms of the OPE these modifications amount to the presence of 
dimension two operator in (j5J); this means that confinement effects manifest themselves 
already at relatively short distances. 

A simple way of realizing such an approach is to introduce a small tachyonic mass 
into the Coulomb gluon propagator [13]: 



= ~^—-2 ^-2+~V ( 3 ) 

—A 2 + q 2 q 2 q 4 

where A 2 > is a constant; we assume q 2 ^> A 2 ~ Aqcd- It is easy to see that a 
Fourier transform of this propagator V(r) ~ / d 3 q exp(iqr) G(q 2 ) yields the potential 
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of the form (pQ). While this may seem a completely arbitrary procedure, the number 
of options one can utilize to modify the analytical structure of the gluon propagator is 
quite limited: indeed, confinement means that the gluon propagator cannot have a pole 
at q 2 = 0, or any spectral strength at q 2 > in general; any modification of spectral 
density at large — q 2 3> Aq CD would spoil a well-established perturbative behavior; and 
any structure away from the real axis would contradict unitarity. We are thus left only 
with a possibility of spectral density supported at small and negative q 2 < 0, and the 
pole of (EJ) is a simple way of achieving this. 



3. J/ip production at small x: strong color field effects 




Figure 1. One of the diagrams describing J/ip 
production in pA collisions at high energy; from 



Heavy quarkonium is characterized by three dimensionful scales: heavy quark 
mass mh, bound state radius R and binding energy e# (in perturbation theory, 
R ~ (m^as)^ 1 , tB ~ niha 2 ), which satisfy the hierarchy ^> R^ 1 3> e#. 

The size of light hadrons is ~ Aq CD , 
and the strength of color fields in- 
side them is E 2 ~ A QCD . Be- 
cause mh ^> Aqcd, the produc- 
tion cross section of heavy quark- 
antiquark pair can be written down 
in a factorized form as a convolu- 
tion of partonic sub-process cross 
section and the parton distribu- 
tions inside the light hadrons. The 
corrections to such a factorized 
form are suppressed by powers of 
(Aqcd/^/i)™ <C 1. The situation is much more complicated for the case of quarko- 
nium production, because it happens over time scale of t/ orm ~ e^ 1 which is still short, 
tform < Aqcd but leaves the possibility of significant power corrections - indeed, the 
description of J/ip production in hadron collisions has been notoriously difficult. 

These effects are even more pronounced in nuclear collisions and/or at small Bjorken 
x, when the strength of the color field increases significantly and is characterized by the 
saturation scale Q 2 (x\A) ^> Aq CD increasing with the nuclear mass number A and 
1/x: Q 2 (x;A) ~ A 1 ^ x -0 ' 3 . In this case the powers of {Q s /tB) n and (Q s R) m cannot 
be neglected and have to be resummed. Diagrammatically, these powers correspond 
to diagrams such as the one showed in Fig. 1 [14] . These diagrams describe the 
broadening of heavy quark momentum distribution due to multiple scattering inside 
the target. Because the momentum distribution inside quarkonium is quite narrow 
(with relative momentum ~ R^ 1 ) such broadening makes it more difficult for heavy 
quarks to bind, reducing the quarkonium production cross section - the mechanism 
considered earlier in Ref. [15J. The universal behavior of parton densities at small x 
encoded in the dependence of saturation scale on x translates into an approximate xp— 
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Figure 2. The approximate "xf scaling" of J/ip production in pA collisions; from [14J. 

scaling of J/ip production in pA collisions |14J . The increasing shadowing at small X2 
(nuclear parton's x) is a generic feature of all approaches to nuclear modification of hard 
processes; however the power corrections to the J / ip production break down the scaling 
in x 2 - The factorization of the nuclear modification factor into a part coming from the 
nuclear parton distribution and a part describing the final state absorption is severely 
broken at high energies and/or forward rapidities. 

4. J/tp in the quark-gluon plasma 

Perhaps the most striking theoretical development of recent years is the lattice 
observation of J/ip survival in quark-gluon plasma at temperatures well above the 
deconfinement transition temperature T c , up to about 2T C [TBI [T7] . This result is in 
stark contradiction not only with expectations based on perturbative expressions for 
the screening Debye mass, but also with the calculations of quarkonium dissociation 
rates by thermal gluons [T8J [19] (quarkonium break-up in this case is analogous to 
photo-effect in atoms [20J; this scenario should be relevant at strong coupling, when 
the binding energy of J/tp in the plasma is large compared to the temperature |21j). 
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Figure 3. The behavior of gluon condensate as a 
function of temperature; from 



The basis of these calculations is 
again the OPE: the dissociation is 
described by the imaginary part of 
the quarkonium scattering ampli- 
tude [22], whose real part is re- 
sponsible for the energy shift of 
quarkonium energy levels as de- 
scribed in the vacuum case by the 
quadratic piece of the potential 
(J2J). In the case of quarkonium in 
the plasma, the leading operator 
which is responsible for the mass 
shift - the gluon condensate - be- 
comes temperature-dependent. In 
hadronic phase, one can evaluate 
this temperature dependence by us- 



ing the relation 
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Y,9h j 



d 3 p 



-n(E h ){h(p)\-G'\h(p)), 

1 7T 



(4) 



where n(Eh) are the thermal distributions of hadrons h with degeneracy factors gh', the 
matrix elements of the scalar gluon operator over hadron states can be evaluated with 
the help of scale anomaly. The result |23j is shown in Fig. 3: even though one can see 
that pions and even nucleons do not affect the condensate significantly below T c , the 
presence of massive hadron resonances as described by the Hagedorn spectrum causes a 
sharp reduction of the condensate in the vicinity of the phase transition. Such a change 
would certainly affect the mass of J ftp - why is it not seen in the lattice data? Moreover, 
lattice results indicate a rapid variation of e — 3p above T c ; since this quantity is related 
to the gluon condensate, one would expect a change in J / ip mass also in the deconfined 
phase - again, a prediction not borne out by the lattice data. What went wrong? 

The puzzle can be solved if the dimension two operator discussed above plays a 
dominant role in quarkonium dynamics, and if its expectation value does not change 
significantly at T < 2T C . In physical terms, this would mean that residual "short strings" 
survive well into the deconfined phase. To get a rough idea of how this scenario might 
work, consider the Coulomb propagator ([3]) modified by the finite-temperature screening 
effects [23] : 

where the component n o of the polarization tensor reduces in the static limit to 
Debye mass squared. The Fourier transform of (jSJ) yields a potential which exhibits 
a residual confining interaction even above T c , see Fig. 4a). One can define also an 
effective coupling [21] through a(r,T) = |r 2 dV ^ T ^ ; the result is shown in Fig. 4b). 
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A comparison of the results of Fig.4 to the corresponding lattice results [24] reveals 
a qualitative agreement, surprising in spite of the crude model we used. Recently, 
thermodynamical properties of the plasma have been investigated in an approach with 
"Coulomb confinement" in the gluon propagator, with encouraging results [25] . 

It is also well known that retardation effects play an important role in quarkonium 
dynamics; a sharp way to elucidate them is to compare correlation functions evaluated 
directly on the lattice to the ones computed by using the lattice potential as an input [26] 
(see also [271 EH EH] ) • Since confinement effects as modeled by the Coulomb propagator 
(J3]) are instantaneous and the screening effects are not, we expect that at short time 
scales the response of quark-gluon plasma will be similar to a response of a confined 
system even if at large time scales the color charges are deconfined. One may speculate 
that this effect is at the origin of both J /if) survival in the plasma and a strong quenching 
exhibited by light and heavy partons. Regarding J /if}, in this scenario it survives in the 
plasma despite being bombarded by thermal gluons because the heavy quarks are still 
bound by a confining potential at temperatures T < 2 T c . 




All of the discussion above assumed that the quarkonium is at rest; how would the 
picture change for a quarkonium traversing the plasma at finite momentum? At weak 
coupling, the answer has been known for quite some time [301 ED] : the gross effect is the 
decrease of the screening length caused by a larger parton density in the rest frame of 
heavy quarkonium. Similar result has been obtained recently for M = 4 SUSY Yang- 
Mills theory [32] which can be solved even at strong coupling: the screening length is 
reduced by 1/a/7, where 7 is the Lorentz factor of the moving pair. Experimentally, 
a stronger screening at finite momentum has to be disentangled from the suppression 
caused by the gluon fragmentation contribution to J / if> production expected to dominate 
at pt > 3 -T- 5 GeV - since the gluons are known to be quenched, so should the J/^'s. 
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5. Why are RHIC and SPS results alike? 

High statistics data on J/ip production in nuclear collisions are available at 
present from PHENIX [331 M, GS] and NA60, NA50 Collaborations [361 EZ] 
at RHIC and SPS energies. A surprizing feature of the data is the simi- 
larity of nuclear suppression seen at mid-rapidity at RHIC and SPS in spite 
of a large difference in energy densities achieved at these vastly different en- 
ergies. A higher density of gluons produced at RHIC was expected to 

result in a stronger suppression 
- but what if, as we discussed 
above, J/ip survives well into the 
deconfined phase? The survival of 
J/tp up to ~ 2 T c seen in the lattice 
data means that it can survive up 
to the energy density e ~ 16 e c 
significantly above the critical one, 
e c ~ 1 GeV/fm 3 . Given that the 
initial energy density in central Au- 
Au collisions at RHIC is estimated 
at e inU iai - 20 GeV/fm 3 (see e.g. 
[38l [39]). it is conceivable that 
even at RHIC the directly produced 
J/^'s are still largely intact, and 
Figure 5. J ftp suppression at RHIC and SPS; the observed suppression stems 
from [30] ■ from the depletion of ~ 40% of 

J/^'s originating from the decays 
of excited x an d states. Such a 
scenario was put forward in [30]; it seems to result in a reasonable description of the 
data, as seen in Fig.5. An alternative explanation invokes the regeneration of J/i/j from 
charm quark pairs [3H 1321 133] . Since the recombination is quadratic in the density of 
charm quarks, the signature of this mechanism is a narrowing of rapidity and transverse 
momentum distributions of J/^'s produced in nuclear collisions. As we discussed above, 
a narrowing of rapidity distribution in nuclear collisions compared to pp ones is also 
expected due to shadowing effects, which suppress J j ip yields at forward and backward 
rapidities. We now need high statistics dA data to evaluate the contribution of this 
effect to the observed rapidity dependence of J/ip suppression in AuAu collisions at 
RHIC. Another important measurement is the azimuthal anisotropy of produced J/ip's. 

I am indebted to Yuri Dokshitzer, Frithjof Karsch, Eugene Levin, Agnes Mocsy, 
Marzia Nardi, Helmut Satz, Kirill Tuchin and Valentin Zakharov for enjoyable 
collaborations and stimulating discussions on the topics presented in this talk. This 
work was supported by the U.S. Department of Energy under contract No. DE-AC02- 
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